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Summary 


The differential equations of quantum electrodynamics and meson theory are formally 
integrated in the higher approximations without use of the interaction representation. It is 
first shown, that the commutator formulae of ScHwrncreR and the P-symbols of Dyson 
satisfy the recurrence formulae in the HEISENBERG representation, and the explicit calcula- 
tions needed are at least not more complicated than before. This method of integration, 
however, has the disadvantage of not avoiding the complications involved in the normals of 
the space-like surfaces, as we have to make an explicit use of the interaction Hamiltonian. 
The real advantage of the HEISENBERG representation is thus lost. In the second part 
the differential equations are therefore solved directly by being substituted into themselves. 
Explicit formulae are given in the simplest cases and the analysis is carried so far, that the 
general structure also of the more complicated expressions can be inferred. The calculations 
are at this point a little more delicate than the use of the Sy (x) and Dp (x) functions (at 
least for particles of spin $, where no surface terms appear) but give, on the other hand, an 
interesting insight into some of the prescriptions of the earlier method. At last the theory 
for particles of integer spin is also investigated, and the similarity to the case of spin 4 is 
stressed. In the appendix the removal of a special form of indeterminacy in the theory is 
illustrated by a simple example. 


1. Introduction 


In a previous paper! the author has developed a method of solving the 
differential equations of quantum electrodynamics (including the renormalization 
of charge and mass) without using the interaction representation. In (I) a few 
simple examples were discussed, and it was shown that the formulae obtained 
agreed with those derived by ScHwiNncEeR and others. In this paper we want to 
examine the formal solutions a little more closely in the higher approximations 
and to compare the results of the investigation with the calculations of 


1 G. KAriin, Arkiv for fysik Bd. 2, Nr. 19 (1950); here quoted as (I). 
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Dyson.”)? It will appear, that our solution is always real, and that it agrees 
with the formulae of Dyson, when no real processes can take place. 
The method referred to above consists essentially in expanding every operator 
(in the HEISENBERG representation) in a power series in the coupling constant 
and to substitute the expressions into the differential equations. If the coef- 
ficients of every power of e is then put equal to zero, we obtain ((I) equations 

(25) and (26)) 


a ay y n+1 
(v5 as m| yer) (x) = ; 2 {A? (2), yy yr—™ (x)} ae Pass yinti—-™ (9) (i | 
DAr @) = — 5 DY), x. @)I- (2) 


These equations are the starting point for the discussion given below, and we 
want to study their formal solution in some details. 

The notations used are essentially the same as those of SCHWINGER and Dyson. 
We thus denote the space-coordinates by (x1, %2, 73) and use %4=1%) =7¢ as time 
coordinate. The four-vectors will not be distinguished in any special way, as 
in any case they are generally easily recognized. The scalar product between 
two such vectors a-and b is simply written as 


ab = ay by == dg be ae az bg — Ap bo. 


However, if we want especially to emphasize the spatial part of some four- 
vector, we denote it by a bar 


@ = (ay, dz, dg) and ab = Gb — ay bo. 


The four-dimensional integration element is, as long as no error may arise, 
simply denoted by dz. An expression of the form 


e J Epa 
| i p2 + m 
is thus to be interpreted as 


fff fe F F 1 e% (121+ P22 + Ps X3— Do Xo) 
a aed Be Ube 

: POPU DMOO ES get we pat me 
and so on. 


As the diagonal elements of the metric tensor are all equal to one (and the 
non-diagonal equal to zero) we do not distinguish between covariant and contra- 
variant indices. 

The matrices y, are the well-known ones, with four rows and four columns, — 
fulfilling 


> F. J. Dyson, Phys. Rev. 75, 486 (1949). 
* F. J. Dyson, Phys. Rev. 75, 1736 (1949). 


372 


ARKIV FOR FYSIK. Bd 2 nr 37 
Yo Yu + Yuyy = 2Quy = 20uy 
and e.g. the expression 


yp (x) 


has to be read as 


Dees a ai 


etc. The A, (x) operators are the electromagnetic potentials, the w(x) operator 
the usual Drrac spinor and y* (x) its Hermitian conjugate. The quantity yp (2) 
is defined as 


The commutator 


stands for 


y (x) Vr? (x) ry (x) ve y (x) ca és (Yr)aB (Pu (2) pe (a)— peg (x) Wa (a)) 6 


a, B 


The other notations are either obvious or explained in the text. Natural units 
(i =c = 1) are used throughout. 


2. The singular functions 


In what follows we will frequently make use of the various singular func- 
tions defined by ScHWINGER and others. For the convenience of the reader a 
short exposé of these functions will be given here, although we have nothing 
essentially new to say at this point. 

We begin with the function A (x) defined, more or less arbitrarily, as the 
solution of 


a , 1 
(= nt)Z (@) = — 32) = — By | evtap (3) 
which is given by 
7 ) re p {a dip 4 
A (w i< (2s)4 p+ m ( ) 


(compare ScuwincER*, appendix). The letter P in front of the integral sign 
indicates, that the principal value of the integral has to be taken in the points, 
where the denominator vanishes. The p-integrations have to be performed in 
the following order. The p,-integration is done at first, and the principal value 
is taken at the poles. After that, the integrations over the three-dimensional 


4 J. ScHwincER, Phys. Rev. 75, 651 (1949). 
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angles are carried through and, finally, we may then integrate over the absolute _ 


value of the three-dimensional vector /. 
It is now convenient to write 


A (2) = $[Ar(z) + Aa(z)] (5) 


where the two new functions Ar(z) and Ax4(z) are given by 


1 om 
Cr 
1 4p 
ae peut aoe ioe (7) 
C4 


The  -integrations are here performed along the two complex contours shown 
in figure 1. 


Cr 
| 4 
= prem? - pit m* 
CG, 
Fig. 1. 
We will also make use of the functions 
aN a (22)* J p2 + m? ~ (8) 
# 
and 
ibe! give 
Oye w= _ ff 
AN @) a (22)* lee + 7a (9) 


where the paths C and (©, are illustrated in figures 2 and 3. 
The factor ; in (9) is inserted to secure agreement with the notations of 
SCHWINGER. Obvious: all these functions are not independent, and we have 
Ar(x) = A(z) — $A(2) (10) 
Aa(t) = A(@) + $A(@). (11) 
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Fig. 2 Fig. 3 
Further 
(4) — m*) An(z) = (C1 — m*) Aas) = — 6 (@) (12) 
Se es (lt) AO (2), (13) 


Due to the factor e** = eikz-ikot the integral over a circle of very large 
radius will disappear in the upper half-plane if x) <0 and in the lower one 
if z > 0. This means that 


Aa (a) = 0 (14 a) 
A («) = —$ A(z) (14 b) 
Ar (2) = —JA (e) (14 c) 
if r > 0 and 
Ar(z) = 9 (15 a) 
A (a) = 4 A(z) (15 b) 
Aa (x) = A(z) (15 c) 
if 7) <0. Hence 
A (2) = — $e (a) A (2) (16) 
where 
Xo 
é (zx) =e aol (16 a) 


From the integral representations (4) and (8) we further infer, not only that 


A 2) = — A(z) (17) 
A(— 2) = A) (17 a) 

but also that 
A(#, —%) = — A(#, %). (18) 
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Hence 
A(z, 0) = 0 (19) 


which, due to reasons of invariance, can be generalized to 
A (a) = A («) = 0, when a? > 0. (20) 
Formula (20) is necessary to make equation (16) truly invariant.* 
The integrals (8) and (9) acquire contributions only from the poles inside the 


contours and, as a simple investigation shows, these functions can equally well 
be represented by the formulae of ScHwINGER.* 


G 
Fig. 4. 
= ae ip x 2 2 ¢ 
A(t) = Gaps { ave 9a" + m*) (0) (21) 
| : 
(1) ee ee tp x 2 ac 
DG) @ wpe | ave 6 (p* + m?) (22) 
As last we also introduce the function 
at etpt J 
Ar (a) = i(2 a)* | p+ mm P (2a 
Cr 


which has been extensively used by Feynman? and Dyson. 3 
The functions S(x), S(x),... etc. are given by 


(2) = (yg, —m) Ac) 


° R. F. Feynman, Phys. Rev. 76, 749 (1949), ibid. 76, 769 (1949). 
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Se er Oe fee a ee es 


and D(x), D(z),... are obtained from A(x), A(a),... if we put m = 0. 


3. The solution of Schwinger 


In this section we want to verify that the formal solution of ScHWINGER 
(compare e.g. Dyson? equation (25)) satisfies our recursion formulae (1) and (2). 
Although this is more or less evident, due to the fact that the two expressions 
have been derived from the same starting-point, it may be of some interest to 
consider the detailed calculations necessary. For this purpose we require a 
theorem, which can be stated as follows. 

Ii A(x), B(x) and H(z) are three operators, and if we define two sets of 
new operators according to the scheme 


5 1 
A™ (x) = i* i dz’ i) (Nr ee f Cre Ni (ee) eee [ae a) ee 
: and 
aa} 


BO (x y= far Peer feo... HO), BOT (25) 


we have identically 


yn 


L(A@B@)™=% fae yee Cet (2), (oe ala (oy, 


—o 
n 


A (x) B(a)]... J) = >, A™ (x) Bo-™ (x). (26) 


m=0 


To prove this formula we expand the iterated commutator in the following way 


[H (a), [...[H (2’), A (x) B(a)].. J] = [2 (@"), [-. [2 (@’), A(@)].-. 


. B(x + Daw A (a**}), [H (a'-), [.. . [H (v’), A (a)]... 
. J). - LA (@’), B(@)] + Pies [... [H (e**), [7 (@**), [.-. 
. [ (@i+), (A (ai-9), [... [A (’), A (a)]... TT. 
-(H (x'), [H (x), B (x)]] + ae ee (27) 
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The first term in (27) gives, when integrated, directly A” (x): B(x). The second 
term gives a contribution of the following form 


Xo X' Nite n | 
(eae ] a2" .7 f da” BD LG ba ae cae © ee (28) 
= Sey — 0 i=1 


This expression, however, is equal to 


gn—-2 


[aw ei eu i be OS (eta eae faa g (0) (29) } 


which is easily shown from n to n+ 1. Supposing (28) and (29) to be equal | 
up to some value n, we get 


n+1 


Xo X'o xy 
[de’ [ dz” oy i} dg > SO a a ree) = 


ee anal 
= faan if dese i da” f(a"... 2") g (atl) + 
Xo Zio gn-1 eg 
+ fda’ [ NOE Sa i) dx" f(x”... x) [dart g (or) = 
Xo aon eo 
=| da’... if dep Gh 2) | degen), (29 a) 


The first term in (29a) is obtained from 7 = 1 only through a renumeration 
of the variables. The second will appear, if we apply our theorem (which was 
supposed valid up to the integer m) for the variables x’’,...a"*!. As (28) and 
(29) are equal eg. for m = 2, they are thus equal for all integers n. 

In a quite analogous way we obtain, more generally 


ie gnr-1 
ro ) 
[da f dar >; f(a” gutlgh-1 girth tpl a") 
aoa -—« i> tg> +++ >%, 
az) uty rah a Xo 
“gat ee) | das. | eas GRE aay ae 
es ao soe 
qgn-1 
vO 
: [ da" g (ee phe gents), (30) 
—o 
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Equation (30) is proved by observing that if it is valid for some integer r, 
we have 


gnr—-1 
Xo oO . 
{ da’ erat | da” = bE eed 8. ar er are gtr ttl Get a). 
— oo —oo > se >the 
n—-?r Xo air+i—t etrt+1 
-g(a*... art1) = »} fda’ 4 F datr+1 | dafrt+iti... 
tp 41=1 —00 08 —00 
yo art 
[ Ce te Go test art te) if Cie Lae 
—oO 
gr—1 
0 
f age ht oto tere). (31) 
—o 


The right-hand side of equation (31), however, is of the same form as (28), 
and we infer with the aid of (29) that (30) is valid also for r+1. Forr=1 
formula (30) is proved earlier, and it is thus correct for all integers r. If we now 
integrate both sides of equation (27) over the 4m-dimensional “‘tetrahedron’”’ 


Xo x'g ay? 
[ da’ fda’ f da” 
—oo —o —0 


we obtain from each term of the right-hand side an integral of the form (30), 
that is, a product of two operators A™ (xz) and B’-™ (x). Obviously, all values 
of m less than n will appear, and we get 


n 


(A (a) B (a ) 7) 2 A™ (x ee wig ) 


as desired. 
After this preliminary calculation we have no difficulty in showing, that the 
expressions 


yO (2) = fda’... f da (He), 6... [A (@) ye]... (82) 
and 
ay gat 
a | ete ee (a) toe Ha), A, (z)\. (33) 


satisfy the recursion formulae (1) and (2), if the part containing the mass- 
renormalization term is dropped, and if H(z) is given by 


H (x) = — 3, (a) Ay” (2). (34) 
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We first consider the equation for y(x) and obtain from 
x 


[ure y@de = — fae a? w) S119 @), oy @) v @)] = 


—-3o 


= [da S(a—2’')yA® (x’') y© (a) da’ = 
+00 
=— f Sr(e—2')y A" (wv) p® (2) da’ (35) 


immediately 
gt 


(, ee + m| yr) (x) = ma [av o(e—e) | Che cae , {dom LH (o™), 


[LE AP ov) = 5 DAP vA}. 86) 


m=0 


The calculations for A, (x) are similar. We have 


J da’ [H(2’), A (2)] = -if [ip (0').y p® («’))- 
-[A,” (a'), dp (2) da" = — | [Pen e.n9 © (eda = 
1 +00 : 
= 5 | Pele — 2) f° '), x29 de (37) 
and hence 
nm+1 vs ras 
mae ian = | d(x — w') da’ | dx’ 


J dat) CH (a), [.. [1 (2), (P (2’), uy (w)I) . = 


ge pI [vy (a), yu pe” (a). (38 


The equations (32), (33) and (26) together supply the complete formalism of 
SCHWINGER. 


380 


ARKIV FOR FYSIK. Bd 2 nr 37 


4. The S-matrix of Dyson 
By starting from our equations (1) and (2), we want in this paragraph to 
obtain the solution given by Dyson.2 For this purpose it is most convenient 


to calculate the elements of the S-matrix directly. This quantity satisfies the 
following equations, given by Yana.® 


[S,p (z)] = — eS [{ S(x —2')yA(z')p(a’) da’ (39) 


[5,49 @y] = ‘fs 58 [De Diz 24 (@), av @nae. (40) 


For the derivation of these formulae the reader must be referred to a forth- 
coming paper by YanG and FrtpmMan. Expanding all operators in a power 
series in e, we have 


[S@+D, y (2)] = —4 D S™ [8 (e—2') (A, @)yp@yr-™ da’ — (41) 
m=0 


[S@+D, A! (2)] = ; sm [De (2 —2')[y(e'), ya p(w irda’. (42) 
0 


m= 


Dyson’s solution is 
Sm OF f fdas. dae P HG)... He’) (43) 


where H(z) is given by (34) and the P-symbol is defined in the following way” 


{ A (x) B(x’) when x, > % 


) B(x’) = 
P(A @ Bix’) | B(x’) A (a) when 2 > 2,. 


(44) 


The generalization for more than two factors is obvious, From the definition 
of the P-symbol we find, without difficulty, for arbitrary operators H (x) and 
F (2) 


fe [dadz' Se Oe A) oe te) oa) = { Nee [dada.,.dz: 


mee % 7 . 
-P(H(a")... H(e’)) F(a) —n f da f da'{ fda"... fda" P(H@ 


HG’): (H(@), F(a —(n—1) [ do” [ fan”. . { dat P(H@"), 


§ Compare (I) note 3 and a forthcoming paper by YANG and FELpMAN. 
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o!) [H (a), [H (2'), F (a) —-- I} = 
eft der. da PH Ge) HOO) keh 


—nf de esate ... [ da®- P(HG@")... H@")(H a’), F(@))+ 


n(n — 1) { de iho faa" fae” 7 fda 


~P(H(a")... H(@')- [H (#”), [H (@’), F(a) — + a | 


If we substitute equation (43) into the left-hand side of (41) and apply (45) with 
F (x) = y A® (x) y® (x) we obtain, using the results of the previous section 


eee |e le 
aa [- faz .. dat} (n +1) P (Ha)... H@)[H@), y® @) = 
(n 


oy Ne 
peed fo fae... dams ea) PG)... Hla A ay @) = 


n! 


= —1 >) S™ [dx S(e—a')y, (Aa) y@yr—™. (46) 
m=0 


Substituting (43) into (42), we obtain in the same way 


(= Oa n+1 n+1 wr peel as (0) Foam = 
Ani da’. . de OPW a)... HG) [A @), A, ain = 


a, 


4) fda, Ba is Dax P(A Ge) .. Ae’) [p® (x") 


Yew) = 5D) SM fda! Die — a) 9), Ya wl. (47) 


The calculations given in the last two sections verify that the integration of our 
equations (1) and (2) is exactly equivalent to the solution of the equation of 
motion in the interaction representation’ 


,0 o|> 

lvlsl> — w(e)|y[o>. (48) 
(The notations are explained in the work by Scuwincer.’) Our problem can, 
of course, be considered as solved, when we have obtained the expressions (32), 
(33), (26) and (43). In practical calculations, however, the iterated commutators 
give rise to rather lengthy computations, especially concerning the ¢-functions. 
The P-symbols of Dyson are much more simple to handle and the ¢-functions 


” J. Scuwincer, Phys. Rev. 74, 1439 (1948). 
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are there cause of considerably less trouble. Both methods, however, make use 
of the interaction Hamiltonian H (x), and although this fact is no serious draw- 
back as long as we are only considering particles of spin 4, a quite different 
situation arises, when particles of integer spin are also treated. In this 
case the Hamiltonian contains terms, which are dependent of the surface o (in 
the notations of SchwineER) and thus not truly invariant. As a matter of 
fact, these non-invariant terms will all drop out of the end formula but only 
after a long and comparatively complicated calculation.2 The disappearance of 
the non-invariant terms is closely connected with the derivatives of the e- 
functions, the occurrence of which in the formalisms of ScHwincER and Dyson 
comes from the boundaries in the z-integrations. The integrals of equations 
(24) and (43) are effectively invariant only if the operators H (x) commute on 
every space-like surface, and this is not the case if the interaction energy con- 
tains derivatives of the field operators. The mission of the non-invariant terms 
of H (x) is only to counterbalance these inconvenient properties of the integrals. 
On the other hand, it has been shown in the paper referred to earlier! that, 
for the case of scalar particles in interaction with an electromagnetic field, it 
is possible to make calculations in the HrEISENBERG representation and obtain 
invariant results directly without reference to any system of space-like surfaces. 
We therefore wish to solve our equations in a manner more similar to the 
procedure used in (I) and do not want to make use of the integral formulae 
above. 


5. Direct integration of the differential equations of quantum 
electrodynamics 


As has been proved earlier!, we can construct a solution to our equations 
in the following way 


ytd) (x) se = | ic (x —z’) ne? {4™ (x’), Vy ar—m) (x')} da’ ies 


© 
m=0 


m=2 


n+l 
— [se (a — 2") >) dm™ yrtt—m (¢')da" (49) 


AN?) = 5 | Dalo—2!) due Dy (e’), roy (eV dx". (60) 


m=0 


Equations (49) and (50) permit us to compute the n + 1*t approximations of 
p(x) and A, (a), when the operators A\ (x), p® (a), ... AY (x), y™ (x) are 
already known. The constants dm” have to be determined in such a way 
that all remaining divergences can be interpreted as charge-renormalization 
factors. The mass-renormalization is thus completely absorbed in the (infinite) 
quantity Om. As is easily shown! 6m?"*t” always vanish, and the summa- 


8 Pp. T. MarruEws, Phys. Rev. 76, 1657 (1949), ibid. 76, 684 (1949). Compare also R. 
Jost—J. Raysxi, Helv. Phys. Acta 22, 457 (1949) and G. KAtimn, Helv. Phys. Acta 22, 
637 (1949). 
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tion in the last integral of (49) is effectively restricted to even integers m. | 
If m is even, m takes the values 2, 4,...m, and if m is odd, the corresponding — 
set will be 2,4,...%+ 1. The solutions given above fulfil, from a formal point — 


of view, the following boundary conditions 


lim y™ (a) > 0) a% > — oo 
51 
a (1) 


lim AW” (2) +0) n +0. 


The incoming fields (a) -> — co) are thus described entirely by the terms y (x) || 
and Aj’(x). It is, however, also possible to construct other solutions to the | 


differential equations, using either the advanced functions Sy () and D4 (x) or 
the expressions S(xz) and D(x) instead of the retarded functions Spr(x) and 


Dr(xz). In those cases the operators y® (x) and A‘(x) will have another | 
physical meaning. If the advanced functions are used, they will describe the | 


outgoing fields '(z > +4- co) instead of the incoming ones, and if the solutions 
are built up with the aid of the functions S(z) and D(x) the corresponding 
operators are given by half the sum of the incoming and the outgoing fields. 
As has been pointed out by Yane®, both the incoming and the outgoing field 
operators will fulfil the usual commutation relations (equations (57) and (58) 
below) and are thus connected through a canonical transformation® 


aplout) (a) ea S-ty™ (a) S 
AS" (g) = S-1 AM (2) S. 
The matrix S above is identical with the S-matrix of HersenBera. If we use 


the singular functions S (x) and D(x), we get for the corresponding commu- 
tators 


[4p (2), Ay” (a’)] = 4 ((LAn” (@), 4s” (@')] + LA? (@), AP (@’)] + 
+ [4," (@), 42” (@')] + LAL? (@), Aw). (82) 


As S is unitary, half the swm of the two last terms in (52) is equal to the 
first term, and we obtain the same commutators for Aj)’ (x) as for AY (x) and 
a (x). Hence, we also have 


Ay (t) = Qt An” (2) Q. 


A similar argument gives the corresponding connection between w® (x) and 
wy" (x). From the equations 


Al (2) = AL (a) + Ff Dee — 2) Lp), may @ da’ 


* Compare P. A. M. Dirac, The Principles of Quantum Mechanics, third edition p. 106. 
(Oxford 1947). 
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AS™ (a) = A® (2) + if 2@—-2) pe Py pplede 


we infer, that we must also have 


gor) 14 
wu (z)=Q (x) Q (53 a) 
porn(a) = Ovi” (¢) O 
with the same Q in (53) and (53a). We thus obtain!° 
OF aS: (54) 


If no real processes can take place, the S-matrix and thus also the Q-matrix are 
effectively equal to one, when we want to compute the expectation value of 
some operator. As long as we are only treating self-energy problems it is thus 
possible to simplify the solutions (49) and (50) in the following way 


ye £8 (eee) ) 34 Yr yr™ (a da’ = [8 (g— a’)? 


n+1 


- >) dm™ yrr1—™) (x) da’ (55) 
m=2 
AarD ( y= 5 [Di (ce — 2’) bn [py (22"), yr pr—™ (a’)] da’. (56) 


The operators Aj’? (x) and y® (x) fulfil the following equations! 


[A? (a), AS? (a')] = — du» D(a’ — 2) (57) 
‘pe (x), pp (a’)} = —t8ea(x’ — a) (58) 
< {AP (x), A? (x')} >So = Sus D™ (2 — 2) (59) 
<[we” (a), yp (2')] >0 = Ssu(a’ — 2). (60) 


Equation (59) is actually incomplete as it stands but will always give correct 
results, if the theory is truly gauge-invariant.1' Formulae (57) and (58) have 
already been used several times in the calculations above. 

For what follows it is convenient to introduce the notations 


1 LA te et {Aye 


Qn- if al ny —1 


n— (0) (n\\ 
(© *) Ang (@)} >< fh 


10 Cf, Dyson”, note to section V. 
11 Compare F. J. Dyson, Phys. Rev. 77, 420 (1950). 
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This symbol is actually symmetrical in all the variables x and indices ». This — 
can be inferred, if we consider the difference 4 


{A, {B, C}} — {B, (A, C}} = [¢, [B, 4] (62) 
where A, B and C are arbitrary operators. Putting 
A= AY (o! 
BA 2 


C= A(3,4,...n) 


we get 


Aly 5:3) xen) A Qeds Balmer 
= —16,[4,.+-), D(a’ —2”)] = 0. (63) 


As A(l1, 2) is symmetrical in x and 2’, we conclude from (63) that A (1, 2, 3) 
is completely symmetrical in all the variables x’, x’ and z’’’. From this fact 
the corresponding property of symmetry can be deduced for the iterated anti- 
commutator of four potentials and so on for any value of m in (61). 

The vacuum expectation value of the symbol, defined in (61), is easily com- 
puted according to well-known rules.* The operators have to be taken together 
in pairs in all possible ways, and the expectation value of every pair can then 
be calculated separately according to (59). If m in (61) is odd, the corresponding 
vacuum expectation value is obviously zero. If nm is even, we obtain a sum 
of terms, one of which is given by 


(4)2 6,,», DY (a — 2’) 6.5, D(a” — a") 2.26 D(x"-1— 2"). (64) 


*m—1"n 
The vacuum expectation value of the iterated anticommutator contains (n— 1) 
(n — 3)...5-3 terms of the form (64), one for each way, in which it is pos- 
sible to combine the numbers 1, 2,...n in pairs. 

In our calculations, commutators and anticommutators of the operators (61) 
will frequently appear. These quantities can be calculated successively in the 
following way. We begin by showing that 


4 LAS? (a \ A (2, 3, = maly Las O04 %¥, D ee) Am) (2, 3) eae n) (65) 


m=2 


where 


AM (1,2, 0. ny sal, 2... ed, ee Lee (66) 


| 
Supposing (65) to be valid for some value n, we have I} 


, 


+ £{A(3,... 0 + 1)[Ay (2), Ay (2”)]} = 
-ntl 


a Co) .) 
= ES ase PS ete Lh nye 
D(a 2!) — 5A (3, ...n + 1) banD 2" —2') = 
= — 5 >) dy», D(a — 2’) A™ (2, ... 0 + 1). (67) 
2 


As (65) is certainly correct for n = 1, it is thus proved for all values of n. 
Using this result we are now able to compute 


ee Wee Atl 2 ln) + 
wa4A,. (2), (4, (2), 4 (3, 4, >... n)]] =A 2,2. a) + 
+ (- i) : Orit, D (x — x7) by», D (a — a”) AM (2, 3,... 0). (68) 


+k=3 


Continuing in this way and calculating alternatively the commutators and the 
anticommutators, we obtain the following formulae 
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a (A, (z’), A (2, 3, ...% + 1)] = EA? (2”), [AY (2’), 4 (3, ... 0 + 1))} + 

: 

| mr (1,2; .:.9n), A(m + 1,...n)} = A(1, 2, 2-2") + 

: 


+(- 3) s SoD SAE he EP 


I=j=l FEs= 
=m+1 


"\ 4 m n 
a a 
-AMMOM (1,2...) + (- ;) Pais paler ee Dae 
t<f<k< rs+t+ 
<l=1 = +u=m+1 


- by:9g D (a — 2!) by, », D (at — a*) by,» D(a” — 21) « 


478 
\ 6 
AOD EHO) (6) (W) a2 Vegetal: 69 
A (1,2,3...0) + ( 3) DD (69) 
and 
m4 (1,2,...m)jA(m + 1,274.0 ep > 5y,», D (a? — a) 
4=1 r=m+1 
4 ig m n 
i .) \ ‘ ; 
pa Ls 2, ssi). ot (- :) >; Dd) bx», D(a" —2) by», D (aa) - 
4<j<k=1 r+8s+t= 
=m+1 
\ 5 
v 
- by, », D (at — at) ADMHOOO (1, 2, ny +(- ;) Sy ears (70) 
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The summations over 7,j,... and r,s,... have to be done in such a way, 
that all possible combinations of the differences 2° — x”, where 7 belongs to 


the set 7,7,... and é to 7, 8,..., appear once and only once. 
For the operators yp (x) and y® (xz) the definitions similar to (61) will be 


Ae pee Pi ae Ls eh a RR ae ks rg 


ies = ; we. 
= seri (2'), 81, 2... tn) fp 2%), [9 i, 


Si, 13 te) Gay ee ele ae): 
S (ie 1, ... n) yp (2). JH} (71) 
where | 
Sy 2, 2 nay, 8 @ = a) tn & A pe ee (72m | 


The vacuum expectation value of (71) is easily found to be a sum of (r+ 1)! 
terms, one typical of which can be written 


GA Sp(S® ( — 2) 8 (1, 2,...4) 8" G@— 22) 8G, +1,.-. bs 
2 SP ar — girtl) 8G +1,... 0) Ce 


The other terms are obtained, if the field operators are combined in different 
ways and can contain a product of two or more traces. The numerical value 
of the factor in front of the expression is always 2~’~! (r + 1 is the number 
of S” functions in the expectation value) but the sign will be dependent on 
the number of traces in the product. As a careful inspection shows, the sign 
is given by 


(— tp (74) 


where / is the number of factors in our expression. For the commutators and 
anticommutators between two J-symbols we get formulae similar to (69) and 
(70). Arguing in the same way as before, we obtain e.g. 


etd (1:2, ecm), F (mi 1, Sepa Te te ee eee 
=) (15223 oe me mie Ieee asiain 8) Se Sect ee ee 
at (- ) le Ares acta tae: 
2) 
tas tat Lys. degag 1, 2,58 es ie Se eee le) 
where | 
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Be ty ity Pils Sn a RE Re dy Deeeettn ts rol... 
it tl... n) = © (a’), S(1, 2... ty) {yp (we 
Tyee 5) in Fil? (x’), (1, pee WY) AY (x), 
[--+ = {y® (ate-1), [yp (ate—-141), 8 (ta-1 + 1, ... te) 8 (ie — giatl). 
- S(t. + 1, ... ag) 8 (x's — atp+1) 9 Vga lela ctl ipe (ZF), 
[= [BG Sie +, mp) TU (76) 


Equation (76) corresponds to (68), and the appearance of the two S («)-func- 
tions in (76) is quite analogous to the existence the D(zx)-functions in (68). 
One important difference between the two formulae, however, lies in the fact, 
that the summation in the last term of (75) has to include also the possibilities 
a =. In these cases a trace of the form 


Sp[S (aiet+1 — a’) S (1, 2,... n) 8S (a® — aiet) 8 (ty + 1,2. . ta41)] (77) 


has to be substituted in (76) in the place of the expression containing the 
S(x)-functions. This difference between (75) and (68) has its origin in the fact, 
that the symbols J contain two kinds of operators (y® (x) and y® (a)), while 
the expressions A in (61) are built up from only one type of electromagnetic 
potentials. 

We will make no attempt at writing down explicitly the commutators and 
the other, more general anticommutators between the J-symbols, as they will 
be only rather trivial modifications of the formulae (69) and (70) but will 
nevertheless give rise to comparatively complicated notations. The main struc- 
ture of those expressions is, however, easily inferred from what has already 
been said. 

We are now able to turn to the direct solution of the equations (1) and (2). 
As we do not want to introduce too many complications at the beginning, 
for the moment we drop the mass-renormalization terms and write (55) and 
(56) as 


prt) (x) —= ; | S ( (x aa x’) ) di A™ (& vy yr-™ (a:’)} dz’ (78) 
Ag” (@) = 5 a Dla —a') dE, yay” (eda, (79) 
m=0 


From these equations we get successively 
oe (a) = —if S(e@—z')y A (2’) y® (2) de’ (80 a) 
AP (o) = 5 { Die 2), np ide (80 b) 
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y® (2) = (— iP ff 8 @—@')[8 (1, 2) y (w") A (1, 2) — dyn {py (@’), 


[py (2), Yup (w"’)]} Srv D (a" BB ieee (81a) 


A? (e) = — Pf fT (1, 2)[5un D (e—a') AW (0) + bun, D(a") AW (a')] (81 b) 


y® (o) = (— iP [ff 8 @—e')[5 (A, 2, 3) yp @”) A(L, 2, 3) — 4 8 (1, 2)- 
{py (a th J (3)} (Ao (a’) Oy,9, D (x ” =o”) as Pe (a"’) 8p, 95 D (a’ — a’) ae 
— brn ly® (x), F (2, 3) (Are @") by», Da! — a") + Are (@") Onn 09° 


- Die!’ —a"))\dax' da” da!” (82a) 


A? (a) = ® ff [ {I (L, 2, 3) (bun D@—#') AQ, 3) + dun D@— 2" AA, 3) + 
+ dy», D(a—a'") A (1, 2) —F (1, 25 3) (8y», D @—a") 8y,», D (e"—2!") + 
+ dur, D (@—2"”) 64», D(a'—2'"))} da’ dx! da" + eff D(«—z’)- 
bun Sp [S (x —a'") 8 (3, 2, 1)]6,,», D(a" —a'") da! da’ dx’. (82 b) 


The last term in (82 b) is obtained from the commutator between wy” (x) and 
wy (x). This expression gives 


ating = 
4 [p? (a’), Voy yl (a’)] = CoE f fae dz” {{p (a’’’), § (3, ty 2) y (x’’)]- 
2A (3,254 fp (ey, (3,1; 2) pa} 4, ae (83) 
The first term in (83) is equal to 
— [ {dau da’ J (3, 1, 2) A (3, 2) (84) 
and the second one gives with the aid of (57) and (58) 
tf [da da” Sp[S (a — a”) (3, 1, 2)] br D (2"” — 2"). (85) 
Up till now the formulae obtained have contained only the singular functions 
S(az) and D(a), but in (82b) and (85) also the functions S(x) and D(a) have 
appeared. When we proceed to still higher approximations, more and more 
terms of this form will turn up from the commutators and anticommutators 
n (78) and (79). Computing, as a simple example, the commutator 
4 [py (x) S (1, 2, ... m) A (1, 2, ... mn), 


A(m+1,...n)S(m+1,...n)p (x")] (86) 
we obtain 
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J(1,2,...n)${A(1, 2,...m), A(m+1,...n)} — 
— {Sp[S (a — 2S (1, De NP alt lore ih), 


A(m+1,...n)] =J(l, 2,... n)[A(1,2,...n) + 


+(- ‘) S Dies Dictate) 3y;, D (at — a)- 


t<j=1 r3=— 
=mt+1 


AMMO (1, 2,...n) +---] + Sp[S(a® — 2’) S(1, 2,... n)]- 


|(-3) s S 8y,», D (2 (of — 2) AO, 2,... 0) wel (87) 


t=1 r=m+1 


the rules given above, a sum of terms, each of which contains only S (a), 
S? (x), D(x) and D® (x) functions. The variables of those singular functions 
are always differences between two coordinates «'— a’. If some 7 < m and at 
the same time the attached 7 >m, there is a simple correspondence between 
these terms and the terms obtained from the expectation value of the remaining 
expression in (87). The latter are obtained from the former, if we in every 
term substitute two, four, six etc. “unbarred functions” for the D” and S® 
functions in all possible ways, only remembering the conditions for 7 and ). 
The same correspondence is easily seen to hold even if the more general ex- 


If an expectation value of the first term is calculated, we obtain, according to 
J 
pressions 


Oty Seen) te” (at )\es(e ode. eee.) (88) 


bole 


are used in (86) instead of S(1,...n) yp (x”). 
It is now convenient to split o our operators y” (a) and AY (x) into two parts, 
one of which (Ty (a) and A (@) contains no “unbarred functions”, and 


write 
vy (a2) = Vy (ar) +p () (839 a) 


Aj” (#) = Ay? (2) Ay? (2). (89 b) 


From the definitions given earlier we immediately infer 


eh”) (x) = (— 2) yf ‘fda’. .. da 48 (2— 2 2 2 801, 4): 


ee 


Be Nel Ae oe ee dar eet we nsdn Las 4), (= 1)"s 

° Pah Oven Vary D (x == x) Ove, Ma, D (x + xv") aula fe Ov, oe vet, 

1+ Dy a a 

an (a*2r-1 ae 21) Ala) +++ (427) (1, 2: Aes n). (90) 
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The summation over 7 has to be performed from 0 to 9° if m is even, and 


—= 15 : : ; 
from 0 to — if m is odd. The numbers 7,,... 7%, have to be summed over 


2 
all possible combinations for a fixed value of r. The last sum over the «; has 
to include all ue where each one of the sets 1, ... 243 4, +1,...%5 

.5 } + 1,...m contains at least one of the numbers Cate etme “The 
Gan tions over ie are also restricted through another condition, ‘which can 
be formulated in the following way. Consider one of the sets 74; + 1,... d41. 
To this set belongs at least one value «; and a function D («i — x'k), wherd 
the number a, is included in another set 7% + 1,...%+1. We then say, that 
the two sets 7;+ 1,...%41 and %+1,...%+1 are “connected”. If a third 
set is now connected, in the sense just defined, to one of the two mentioned 
earlier, we say that it is also connected to the other one etc. The remaining 
restriction for the «-summations can then be so formulated, that none of the 
sets 7; + 1,...%41 is to be connected to itself. As an example of a term, 
which does not appear according to this restriction, we can mention 


—4 (a [ff S@— a’) yn fy (e’), J (2,3)} 
A? (2!) bnyeg D (ae! — a!" da’ dao!" dat”. 


(Compare equation (82 a)). 
The corresponding formula for the electromagnetic potentials is given by 


WAL? (a) = (— i f---fada'.. . og 2 Zi Jil? ae. Wachter see 


T ty< 
iinpiget Lon. 1) (ye eS Lene oi) ai a3) 
OY + AD» 


-D (a%27-1 — gf2r) >)’ Ou»; D(a — a) AM) -+-(@2r) (1, 2, ... n). (91) 
j 


Y, v fi 
aH et Re 


The same conventions for the summations over 7, 71, ... 7%, &,... @, are valid 
n (91) and (92). The index 7 in (91) has to be summed over all values from 
0 to n, excepting the numbers a,,... a2. The proof of (90) and (91) is con- 
veniently done from n to » + 1 in the manner used several times earlier and 
will not be given explicitly. In this calculation it is permitted to disregard all 
terms in the right-hand side of (69) and (70) except the first one in (69), as 
those terms will give contributions only to wy” (x) and Ae (aye 

The remaining parts of the operators (?y (a) and ° A (z)) cannot be written 
down in such a compact form, as was possible for ’y (x) and ‘7A™ (a). Their 
general properties, however, can easily be inferred from the remark made after 
equation (87). As these expressions will all appear from commutators of the 
form (86), it is very easy to include, later on, their contributions to the ex- 
pectation values of our operators. We thus turn to the calculation of expecta- 
tion values and choose for this purpose the current operator as a suitable ex- 
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pression. This operator has a fundamental significance as both the interaction 


Hamiltonian and the S-matrix are closely connected with it. If we take only 
the part ‘j(x) defined analogously to (89), we obtain with the definition 


Dj, (x) = >, ent ingle) (x) ; (92) 


the formula 


Ba” (x = (—12) yf: ‘fda’ .. wi Mg > 2 2 | (==1)" 3. 


= Yr . 
ge Lec 0 0, Whom b, .. teh cere. on te Sr ce Ly een)" 
T ge ane T) (a8 en] ee Hp Qoe\ - 
| Ove, 9, 0 (v4 — 2)... Oreerae ov D (x°2r-1 — x27) 
Pine 1D 
pete er (1. 2. 9). (93) 
The summation conventions for r,... 2, are similar to the earlier ones. The 


expectation value is now obtained as a sum of several terms, each of them 
containing functions S"’ (xz) or D® (x). It is at this point most convenient 
to introduce the graphical language of Dyson? ? and Feynman.” If our ex- 
pression contains n + 1 variables z, x’,... a”, we represent it by +1 points 
in a graph, each point labelled with one of the letters a, z’,... a”. For every 
function S(z' — 2’) or S” (a'— a’) in the expectation value we draw a line 
(an electron line) from the point x to the point a’, and for every function 
D (xi — x) or D™ (a — 2x’) we draw a dotted line (a photon line) connecting 
the two points in question. For every operator representing a free particle we 
also draw a line (electron or photon) from the edge of our graph to the point 
that appears as a variable in the operator. In this way we obviously obtain 


a simple diagram representing our expectation value. It is, however, doubtful, 


if this graph will also give a picture of the physical process underlying the 


theory. In the older formulation of quantum electrodynamics the higher radiative 


approximations were described with the aid of virtual photons and virtual pairs 
of electrons. In the papers quoted above, FeynmAN? has developed an intuitive 
method of identifying the electron lines with the virtual pairs and the photon 
lines with the virtual quanta. From the point of view adopted here, this is, 
however, only an extra complication, as we have no need for the conceptions 
of virtual particles, and as the singular functions have the much more natural 
property of being either solutions to the inhomogenous differential equations or 
expectation values. Also in the older formulation of quantum electrodynamics 
the virtual particles were, even from a principal point of view, not observable 
and thus only a mathematical conception. We therefore prefer to regard our 
graphs only as a suitable representation of the mathematical structure of the 
theory. 
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It might be pointed out here, that it is also possible and in some respects 
convenient to represent the field operators themselves with the aid of graphical 
figures. If we use the same conventions as above concerning the S (x — a} 
and D(x'—«’) functions but write a letter py, y or A for every operator 
p (x), py (a) and Aj (2), we can make us a diagram of the “yp (#) and 
©4,,(x) in a simple way. For the operators given by (80), (81) and (82) we 
then obtain figures 5, 6, 7 and 8. The summation rules of equation (90) can 


wrt Agr) 


A 
eo eed f 
7) w x x 
itn Aux) \ 


< 


x rSaeee ee pee al oon es 
Jan x 3 ¥ x 
7 


Fig. 7. Fig. 8. 


here be expressed in the following way. We sum over all possible combinations 
of either n operators A and one operator y or (n — 2) operators A, two oper- 
ators y and one operator y or (n — 4) operators A, three operators y and two 
operators y etc., where the different operators are connected together with 
electron and photon lines. The restrictions on the «;-variables correspond to 
the restriction to only “connected” graphs with no “closed loops”.1? In a similar 
way the ‘'j(«) operators (not their expectation values) can be represented 
by graphs containng y and A as above. In the lowest approximations we 
obtain, corresponding to equation (93), figures 9 and 10. 

We obtain the graphs for the expectation values from the graphs for the 
operators by combining the letters A (or yw and yw) together in pairs and sub- 


* For the definition of these expressions see Dyson® p. 495—496. 


394 


yay 


ao () 
A A 
A 
08 NS, ee ea ee 
§ Yann, Me 
x 
A A 
a A ee ey 
sate j————{ 7 ¥ 
? ¥ ¥ 12 ¥ 
ee oe See 
io ae ~ 
s veae 
Fig. 9. Fig. 10. 


stituting a photon (electron) line for each pair. This procedure has to be per- 
formed in all possible ways, and it is easily seen, that all the connected graphs 
of Dyson™ will appear as a result. No disconnected graph, however, will turn 
up, but these diagrams were shown by Dyson to give only an unimportant 
vacuum fluctuation and no contribution to the expectation value.!? It might. 
rather be considered as being to the advantage of our formalism, that these 
terms will not come out from the calculation. 
For each connected graph of Dyson we will get, not only one term as he 
does, but several terms containing both the S (x), D® (x) and S (zx), D(a) 
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: functions. If we consider a graph belonging to the nt approximation, this 
: 


diagram does necessarily contain exactly ; closed loops, which may or may 


not have one or several lines in common. If we first consider the case, when 
all the loops are entirely separated from each other, we infer from the rules 
given above that the lines joining the loops must be represented by S (a) or 
D(x) functions. This follows from the fact, that the diagrams of the operators 
(e.g. figures 9 and 10) contain no disconnected parts. As a simple example we 


Fig. 11. 


can consider the graph in figure 11, belonging to the one-particle part of 
D7? (x). This diagram contains two loops, one between (z’, x’’) and the other 
between (a, z’’’’). The lines joining the loops are the lines between the points 
(x'’, x) and (x, v'’’). These two last lines are thus represented by the functions 
S(x’’—a) and S(a2—~z’'”). As the diagrams of the operators contain no closed 
loops, one and only one of the lines in every loop in the graph of an ex- 
pectation value must be represented by a S" (x) or D(a) function. It is also 
obvious, that all possible combinations will actually appear, and this means, 
that each loop is represented by a factor of the following general form 
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‘8S (m +1,. 20) 6,5, De = 2). oe 


The modifications necessary, when more or less than one of the lines are photon 
lines, are obvious. For the special graph of figure 11 we thus obtain, apart 
from a numerical factor 
f- ‘fda’ ve Gay (a!) yi (8 (2? — @) D™ (a —2') + 

ae Sh (a G- a’’) D (a —— a’)} 8 (a _— x) viS (x a a’’’) Vo" 

2 {§ (x0! ca ae D” OR ee a’’’) ce gn (a’”" mG a’””") D (2’""" = a’’’)} : 

Yo yp (a’””’) ’ (95) 


Using the following integral representation, given by LutrincEer!® 


6 (a) a 0 (ag) we 0 (An) PMID 


Ag... An a, Ag... An Ay... + An-1 


in —2 


1 @y 
= (= 1 [day | day en. [oda 8°) (ante + 
0 0 0 


+ Ona) nse, — Oot | te Oa Oe La) ay a) (96) 


{the proof of equation (96) is conveniently done from n to n+1) we can, in 
momentum space, write the expression (94) as 


4n—1 
(5 45 | fan Ta O Dat die cite’ —2') +t pea" —a'") + ++ +i Dy _y (1 a), 


I 
245.09. By St) Vans Pog OY eet — tN ae 
1 én—2 
Ons ry (al) [as ate | @ En 31> OY) (he? En_1 + pn—-1 En-2— En) + 
0 0 
tes + pe (Ey — &) + pil — &) + m2 (1 — En_v). (97) 


For the special case of equation (95) we obtain in this way 


lorena ener 
(,",) | far, dhs | ve | Ap... dps e'® PsP») & (py — ky) ya, (ty py — m)- 


"Ya (CY D2 — M)° Yu (Uy D3 — M) Ya, (VY Da — M) Ya, U (Da — he) - 


J. M. Lurrrncer, unpublished letter to the author. 
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1 


SE ta Oe) oe — PH) f apy ate + of + mba —8): 
0 


(p2 + m") (p3 + m?) 


1 


: [ ano @3 n + (ps + m?)(1 — n)) jacil dqdq eit(a'—4). 
0 


[ dkdk' i (9) ya liy@+B—m) yy (iyq—m): yp 


1 


aes pe SAT gets a , 1 U 2, 
(lyq’—m) ya (ty @' +k m) rat (@) ram geawn | Uf (APE + 
0 


+ (q+ hk? + m)(1— 6) { dno (kh ?y+(q' +k’? + m)(1—7n) (98) 
where 
Pe 1 —i@q2Z o,,(0) 
W(0) = Gaye | Ty de. (98 a) 


Introducing the notation 
1 
K (q) = — [dé [dkyliya + k)—m)y.0' (PE+@+h? + m)1—) (99) 
0 


we obtain for the expression (98) 


ry by fasten itya—m_  iyg’— 
(52) ff cada ee ag) Key, 2L—™ Ka u(@’). (100) 


gtm q2+m 


Using the method of Dyson, we obtain instead of the expression (95) the 
formula 


“i 2 
(;) | Mee | dx’ ees aa y (x’) vy Sp (x’ = x") vale (x , x 
, Vu Sp (x as ee’) Ye Spr ae — 7a Vo yo (x pO : (101) 


With the aid of the integral representations of FeyNMAN!4 


2 


=(n—1)! [en fof a Mer 


ay .-- An 


: [an Ln-1 + An-1 ee mF Ln=1) ea dy (Lia) |e (102) 


4 R. P. Feynman, Phys. Rev. 76, 785 (1949). Compare also R. Jost-J. M. Lurrrncer, 
Helv. Phys. Acta 23, 208 (1950). 
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this expression is, in momentum space, conveniently written in the form (100), 
with K (qg) given by 


==> yltiyat+h—m) yn. 
: Jee fetae q+ B+ m8 = 


Introducing the well-known translation of the origin k > k + q(1— &) and ob- | 


serving that 


dkf(k2) (—1)"n 


| 
[ axom (k2 + A) f (K2) = Re Leraye — (104) 


we immediately infer, that the real part of (103) is equal to the expression 
(99). From a comparison of equations (96), (102) and (104) we further conclude, 
that this equality is not limited to the case of only two lines in the loop but 
quite general. If A is never zero in the interval 0 = <1, the right-hand 
side of (104) is real, and the two kernels obtained are identical. However, if 
A does vanish in the interval in question, the integral (n = 2) 


dk 
i (k2 + Aye oo 
Cr 
becomes infinite as!# 
17 1 ee: 
| a te ee ee (106) 
As 
[dk 6™ ( + A) (107) 


gives only the real part of the right-hand side of (104), we have 


1 


eek (108) 


[azom (2 + A) = (—1)*°x(n— 2)! P— 


The vanishing of A will be related to the occurrence of real processes!® and 
gives no trouble in the simple form of diagrams investigated up till now. Both 
the formalism of Dyson and the integration method developed here are only 
applicable, when the elements of the S-matrix are effectively equal to one, 
and it is not astonishing that they will give different results, when this con- 
dition is not fulfilled. 


** Compare Dyson® p. 1744. 
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The method of parametrization used to obtain (100) from (101) with the aid 
of (102) has been found convenient, as it is then especially simple to identify 
the divergences occurring in the theory. It is, however, not impossible to 
apply (102) in other ways too. At first it would possibly seem natural to use 
only one set of auxiliary variables and write for (101) e.g. 


afl: (4) { [ dqdq’ u(q) K (g, q’)u(q’) e*@-9 (109) 


a 
with 


1 Eat UM 
K(a¢’) =5! [de { dry... fdas | [abate piya+ b—m): 
0 0 0 


yaltve—mM) yu(tyg — mM) yo(ty@ + k')—m) yo: 

(gq +.&)* + mm?) 25 + (q'® + m*)(%— 25) + (9? + m?) (xg — a4) + 

+ (q + kh? + m?) (a — xg) + k? (ay — te) + kh’? (1 —a,)}. (110) 
In this formulation, however, the divergences are not so easily separated as 
before. In our formalism, on the other hand, only the parametrization given 


in (98) is natural, and this choice of the auxiliary variables corresponds exactly 
to the parametrization of equations (100) and (103). 


Fig. 12. 


A slightly more complicated situation will arise, when two of the loops in 
the graph have one line in common. As simple examples of such diagrams we 
can use the two graphs in figure 12, both of which belong to the one-particle 
part of j, (x). From the calculations given earlier it can be seen without dif- 
ficulties, that the expectation value of the (x) operator will give a sum 
of terms, each containing one S” (x) or D™ (a) function in every loop. Each 
possible combination will occur, and eg. for the first graph of figure 12 we 
obtain 
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fo fat da 9 pi (8 (0 — 2") 8 (0 —2"") p80" =a) 
-D® (a! — 2’) D® (a'" — a!) + 8 (a! — 2”) yy 8 (a'"’ — a’): 
yy SY (a! — 2’) D(a” — 2") DY (a — 2") + S (a — 2") p: 
SD (2! — 2") y, 8 (2 — 2") D(a" — 2’) DY (a'"” — x’) + 
+ SY (a0 — 2!) py 8 (a! — a!) pS (a — a’) D(a" — 2) - 
-D® (a! — a2!) + 8 (a! — 2") » SY (a’" — 2") y, S® (a — 2’) - 


A ID (a’’” — a) D(a" — a'’) ale S™ (a Mii Me Sigel bs 1 __ gl!) : 
yw S (a —2') D(x" —2') D(x" —2")} 8 (ae — 2) vy, p® (a). (111) 


In these cases we also have contributions from the terms containing the operators 
™y (x) and °'A” (x). The expressions obtained in that way will contain one 
functions S(x) or D(x) in the line common to the two loops and another func- 
tion of that kind in one of the other lines in the loop, not directly connected 
to the point xz Furthermore, those terms will appear with a minus-sign. For 
the graph belonging to equation (111) we get simply 


ee. fe f dz’ be dda yp? (x sa) v8 (a’’”’ ers a’’’) y, 8 (a! = a’) Vy : 
-S (a — 2’) 8 (2' — 2) ppp (2) D(x" — x’) D(x’ — 2"). (112) 


The parametrization of the singular functions in momentum space can now be 
performed in the following way. For the terms containing a S"’ (x) or D™” (a) 
function entirely in the loop directly connected to x, we combine all the ex- 
pressions with this singular function in the same line with the aid of formula 
(96). If the momentum variable in the line in common is called p, and the 
corresponding variables in the loop not connected to xz are p+q', p+q”,...p+Qq", 
we obtain in this way 


1 & n—1 
(—1"fdp [dg fd&... [ d&xd™ (p+ Q24 A) (113) 
0 0 0 
where 
Q = 0" (G1 — fn) seg” (n= aS a (113 a) 
and 
A =m? + Ge (Fn—1 ca En) sh ene ath he (1 = §1) 5 Q? D (113 b) 


The whole expression (113) is multiplied by 


Ba p) te tales 0 (p"’®) aie. 


7112 


ary | ee (pi + m?) p'*(p 


7172 


pp a an") Soe ee) 
6 (p'? + m?) 
TS ae ee mile 2 (114) 
<n (etree) 
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where p’,... p’ are the momentum variables of the loop connected to the point 
xz. We have here made the assumption, that all the lines in the two loops 
whith the exception of two are electron lines. If some more of them are photon 
lines, equations (113 b) and (114) will be slightly changed, but as the modifica- 
tions necessary are quite trivial, we will not investigate this possibility further. 
The remaining terms from Yj‘) (x) can obviously be written 


™n—2 


1 
(—1)" [dp d(p? + m?) [dm ... f dyjn—10 9 (p+ Q'2 + A’) (115) 
0 i) 
where 
O= Ptinza tO (naa — Hn) Fo + (Lm) (115 a) 
and 
A’= q”? Nn -1 a Gs a=) Yel) ane gre (Lin) ma (115 b) 
The terms from ol (x) can analogously be combined in the following way 
1 "nm —2 
(—1)"f dm... f dyn—1f dp d(p? + m2) e(p)- 
it 0 
6") (p + QP? + A’) e(p + Q’) (116) 
with the same notations A’ and QY’ as before. The expression (116) is, if no 


real processes can take place (A’ > 0), equal to 


"pn --2 


a 
OS fame f amar f - 
me Om | OP | e+ erray ON 
i) 0 


and thus (115) and (116) will combine to 


in —2 
2 


im — 1)"-1 | ieee | dma { dy 3 (p? + m?) 0") ((p + Q? + A’) — 


0 0 
1 Serta) | 
4 : Sh (118 
me (p? + m2) ((p + Q’2 + A)" f (118) 
Using the well-known integral representations 
, 1 twa 
6(a) = sy | dee (119) 
p! = es [ aw ey ewe (120) 
OU Divs | w| 


we now obtain 
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%n—2 


Pes Wy W2\ _ =f a om 
(2 We) (13 [w, | ales = dm... dYn-1" 


; fda a—ayrt [dw w1|u| f dpcteurisee a-otmasia'sa naar = 


Nn —2 
= — 4(—1)* p [tea —ar far. a ice 1 f dur wr? er — 
1 at ™n-2 


=(—1)"*-2 [ axa—ay [ an ao [ dm- ot” (B) 
0 0 0 


B= m—Q'? (1 — a)? + (A m2 + Q') (lL — a) = mi? + (ni + - 


+ q!®(1— my) (1a) —Q'2(1— a. 


In the derivation of (121) the transformation 


UW, = Wa 


W, = w(1—«) 


has been used. The new variable-transformation 


"mn —2 
S a fen. pines ding {dp ctmietemestentns area 
0 


(121) 


(121 a) 


(121 b) 


(122) 


will change B to A of (113b) and Q’(1—«) to Q of (113 a). The expression 


(121) can then be written 


's 
Ss 


n-1 


i 
(— taf dé... f d&,8"-(A), 
0 0 


The p-integration in (113) can be performed according to well-known rule } 


and we obtain 
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; —1)"-2(n— 2)! 
| dp d™ ((p + Q? + A) = ( rae ) ae (124) 
The expression (123) is multiplied by 
Mae R LE CMR 0) Sera et Ni (125) 


where p’,...p' have the same meaning as in (114). Equations (114), (124), 
(123) and (125) thus combine to give 


1 fn—-1 


(ee ee ee 
x | ay... f aon love... (ptm) * pt. @H am) 


0 0 


pen-2) (4)| - (126) 


With a simple modification of (96) (some of the numbers a; equal) we can in- 
troduce a new set of auxiliary variables and perform the remaining momentum 
variable integrations without any new difficulties of a fundamental nature ap- 
pearing. The parametrization obtained in this way again corresponds to the 
most convenient one used by Dyson®, and the result of the integration is equal 
to the real part of the formula obtained with the use of the Sr(x) and Dp (z) 
functions. 

At this point it is perhaps not necessary to go further and to write down 
explicitly the expectation values, when the loops have more than one line in 
common, or when more than two loops are mixed together. Also in these cases 
we get one term of the form A~” from the “normal” terms in the loops and 
contributions containing 6”—)(A) and 6”-)(A)-«(Q) from the others. The 


last expression will appear, when more than two loops are mixed together. 
We then get terms of the following form 


[ dp d(p? + m) 8") ((p + QE + AY (ep + Q)— ew) (116') 


Comparing (115), (116) and (121), we conclude 
fdpo (p? + m?) 6—) ((p + Q’? + A’) (1—e(p)e + Q”) = 


1 
= 7 [axa — a)" 6-2) (A (1—a) + ma + Q’2a(1—a)). (121’) 
0 


This result can be interpreted geometrically. To get a contribution from the 
left-hand side of (121), we must have 
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p+ m =0 
(p+ Q) + A’ =0 


Po Pot Qo 


Imo] |p + Qo] 


RP 


(p+Q)?tA'=0 


Fig. 13. 


These equations cannot be fulfilled at the same time, if we do not have (com- 
pare figure 13) 


(Qo|> m+ VA’ 
in all coordinate systems. This can be written in an invariant way as 
Q’?>—(m+VA')P 
which is exactly the condition, that 
A’ (l1—«a) + ma + Q’?20(1—a«) 


will vanish for some value of « in the interval (0, 1). Now, the integral (116’) 
gets contributions from exactly the same values of p as the integral (121’), and 
the only difference between the two expressions lies in the signs. We thus 
obtain 


[dp 6 (p? + m?) d—» ((p + Q’? + A’) (e(p + Q') —e(—p)) = 
= 7% | da (1 — a)"—2  (Q’) 6-2) (A’ 1 — a) + ma 4- Q’2% a (1—a)).  (121”) 
0 
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The restrictions on the graphs introduced above are in no way essential for the 
argument and were only established in order that we should not be burdened 
with too complicated formulae. It seems to be a general property of the higher 
approximations of the present form of quantum theory that, although they can 
be comparatively simply understood from a general point of view, they give 
nevertheless rise to very cumbersome computations, if we try to make an ex- 
plicit calculation for a special problem. 

It now remains to include the mass-renormalization terms of (55) in the 
calculation. This can be done comparatively simply by noting that 6m in the 
equation 


; | 
(> 57 + 1m) vla) =" (Aya), yo (a)} + dm pla) (127) 


plays a role similar to that of iey, A, (xz). Introducing a new vector 


Aj, (2) = Ay (2) — 7-yn- Om (128) 
we have 
0 } j 
(v5 +m) we) = FAs), yo wo) (129) 
0 A, (x) = — "2 oe), yp). (130) 


We thus obtain the solution to (127), if we substitute the total 6m (including 
all powers of e) in one or more of the places for the operators cey, A,” (x) 
in the solution of the equation without mass-renormalization term. These sub- 
stitutions have to be made in all possible ways and the results added together. 
We will not investigate this point further, as the calculation would be merely 
a repetition of the argument given by Dyson. 


6. Particles with integer spin 


In (I) the interaction of particles of spin zero with an external electro- 
magnetic field was studied in some detail. It was, however, also observed, that 
no serious complications will arise, if the electromagnetic field is quantized. In 
this case we have the following system of equations 


+ <4, (x) A, (a), p(a)} (181) 


re f/ 0 . 
0 An(e) = (9 (a) go — °F), w(o)} + 2 Ay (a) (g* Ce) (2). (132) 
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Expanding in powers of e as before, we obtain for the first few approximations — 


(compare (I) p. 190) 


o” (2) = —1 af A‘ (x!) (A. 2") i sh "Ag *)) gp (x')dax’ (133) 


) m 0 0g * (x : 
Ap (2) = —5 it Do—2'){(9* a8; To), y (2)| Bye (134) 
gy? (x) = (—1)? [ [ax dx” E (1, 2) [Ate 2’) 4 a. OA (fc =| ? 


f 
\ 0 x», 
ae 0 0 A (2/— 2" “as 
[ae "xo + At 7 eee ane band ® (a) + 


|e (2’), {(o®* @ — ed |] | (135) 


(2) é , af : ; 0 0 (0) # (7 
a (af / tds i|i[o” (a) Fae ee |. 


[Ree + Ome + 


, 
Tye O ay» 


+ 82") Sunt (2), 6 (2")} | LD ea!) bum A? (2") + 
+ D (x—2"') bury AS” (x')]. (136) 


It is clear, that the general structure of the solutions to (131) and (132) is, 
apart from some differences in signs, very similar to that of the solutions to 
(78) and (79). We have only to introduce expressions of the form: 


[a c= ee si 2a G2) de E (a1 —a") My al 


Vy 0 ay, OE ! Go 
n—1 
0A (e—2’) 
+ t) Pee 
y= [a @ ‘ ) 0 We i O&,, | 
Zz 0 0K (a'-?—2'-}) 
1-2 »t—l1 ae: = a 
[Av aes 3) gt zi OE ee | 


| 
| 
\ 
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y) : a Stihl 
[Rey ae ane |. 


"+2 "+2 


a 0d AR(@"1—a")] . 
[aero 52 + A (a “4 Sent } (137) 


n—1 
| Oa a 


(compare (I) equations (17) and (18)) instead of the factors 


rd C8 We at 8 Cea a oe Pe S (x"-1 — gM) Vn (138) 
and 
0 dp *{a') ot 0 OK (a'—2"") 
(Q) * (in Rage eataits 
tL?" © ane ga | [Ae —2 a5 + aul 
28) 0 0A (a"-1—2") 
ae 2 — 1 (0) (pn 
[Benge +t AG [9 er + 
n—1 
A dg* ©) nahin: 
(0) # , Ci ele Ws 
+S [¢ ) oar eae Meta Pe ae) 
 Sy:r¢41 O(a — 24)... @® ()} seg eric leer (139) 
i<j 
Instead of 
[po (a2), 7nd (e eens (ah 1 — 2") y, p® (2")). (140) 


a special case, written down explicitly in (136)). The analysis of the higher 
_ approximations is now identical for particles of spin $ and spin 0 and will not 
be repeated again. 
It is obvious that an argument similar to that given above can be applied 
to all kinds of mesons in interaction with nucleons or an electromagnetic field. 
If only the operator equations are covariant, we will always get an invariant 
result out of the calculation without any elimination of ‘“‘surface terms’. As 
another simple example we can take the case of pseudoscalar mesons (operator 
p(x), mass m) in vector interaction with nucleons (operator y(x), mass M). The 
differential equations are here 


| {The modifications necessary in (139) for 7 equal to one are obvious and, for 


a a 
(v5 + u) y (x) = a eee, V5» ()} (141) 
fe 
(Om?) 9 (2) =F 9 [¥@): ys ¥(@)] (142) 
Ys = Vi V2yVaya- (143) 
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The recurrence formulae corresponding to (78) and (79) are 


“s (m) (a 
venta) = —F [Se—2) BPTI env ede’ C4 


m=0 


0A (a — a’ s ' ' 
yrrD (a) = ay A= 2) Sy (0'), ¥5 9%» yOr™ (e'y)da’. (145) 


m=0 


It is now evident that we obtain the solutions to (144) and (145) from the 
solutions to (78) and (79) through the following replacements 


0 gp (x x) 
A® (146 a) 
yw (&) > 0x, | 
Yu > V5 Yu (146 b) 
ONO e Ae a | 


D (a' — 2") Oy» > 


, TF. 
Oxy, Oy, 


From equation (146c) we infer, that the expectation values will, in the n™ 
approximation, contain » more momentum variables than the corresponding 
expressions in quantum electrodynamics. The divergences are thus much more 
serious in this kind of meson theory than in the case of particles of spin one 
half in interaction with an electromagnetic field.1® 

The other cases of meson-nucleon interaction can be treated in similar ways. 


7. Appendix 


In this appendix we will investigate a point, which, although it lies a little — 
outside the main discussion, is very conveniently studied in the HEeIsENBERG 
representation. We will not carry through the calculations in all approxima- 
tions, but restrict ourselves to the e?-terms. Let us consider the one-particle 
part of the current operator for particles of spin one half. In the notations 
used in (I) this expression can be written (compare (I) equations (35)—(37)) 


jy (2 -5/ [p (a), YuS (a — 2’) (b a’) — 6 m™ y (x’))], da’ + 


Le, a [(p @!) — d m® p @) 8 (a! — a) Yu, yp (eh da’ — 
“ae | | [py (a’), Ky (a! aye a!) pl’ (da da vac, pol (Age 


The operator ¢ (x) is logarithmically divergent but, using a regularization with | 
respect to the photon mass!’, we obtain 


16 Pp. T. Maruews, Phil. Mag. 41, 185 (1950). 


1” W. Pauvi—F. Viniars, Rev. Mod. Phys. 2/, 434 (1949). Compare also the paper by 
Jost and Lurrincer.* 


408 


ARKIV FOR FYSIK. Bd 2 nr 37 


2 
¢ (x) = ar [ f avax [ae Gyn 2m) ef @+k) (e-2) yy (p') . 


i {d(p2+ m?) 6 (k?) oO (k? + M?’)\ 
U2 (Re + M2) © (p? + om?) (2 + M2) (p? + m2) 


(A. 2) 


This expression can be calculated according to well-known rules, and we obtain 


(x) = $1 (%) + de (2) (A. 3) 


1 uu 
zt M? ; ; u(iyg + m) yp (a) 
= ’ pi a(a—@’) 
$1 (2) aor | auf dof da f ae if q?u(l1—u) + m? (1—u) + Mv GeG 
6 0 


1 U 
at M? nies (2— u) yw (2’) 
SS 4 pig (x—2') 5 
$2 (x) Sor | auf dof aaf ae F gu (1—u)+ m?(1—u)+ Mv —) 
0 0 


As y (x) fulfils the Drrac equation for free particles, we can use 


gtm=0 ° (A. 6) 


in the denominators, and obtain 
if 


1 M? 0 
dy (4) = a2 [udu log (1 + ae =] (rz, m) wy (2) (A. 7) 


0 


1 
m 2 


‘ M 
0 


| 
The infinity (dependence of M) in ¢9(z) can be cancelled, if we define 


1 
m 


om? = 3 | du(2—u) log (1 + ——— 
0 


It is not permitted to put 
( ci ae y= 0 A. 10) 
75a t My ve) = (A. 
in (A. 7) as this expression is multiplied by S (a—a’) and integrated over a’. As 
(y5-, +m) Sea) = — bee) CARLY 


the term 
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[8(e@—2') ¢1(e') da’ (A. i 
is indeterminate and can be given any value of the form 


2 


1 
1 Mu 
— co [udu log (: + aaa) y (a) (A. 1 
0 


with &€ an arbitrary constant. ; CA 
The term containing K,(x’ —x;«—a’’) in (A. 1) also contains a divergency, 
This can be separated in the usual way!’, and we get, after regularization é | 


1 


4 

i : Mu ‘ 

HA = 

° :* 


Collecting the terms from (A. 13), its Hermitian conjugate and (A. 14), we obtain 
1 4 | 
Y 2 

: [ du-w-tog (1 + = sedis ya 2 €) [wy (x), yup (z),- (A. 1 


16 ae m* (1 — u)? 
0 


The second approximation of the current operator thus contains an indeterminate — 
multiple of the zeroth-approximation current. * From a mathematical point of 
view, this corresponds to a solution of the homogeneous equation added to the ~ 
expression given e.g. in equation (49). Such an extra term ought to have been 
avoided according to our boundary conditions (5) and the use of the retarded — 
functions. This is, however, not the case and is probably connected with the 
fact that the operators y® (x) do not vanish for %—-—oco. Hence, the solu- 
tions (49) and (50) do not strictly fulfil (51). According to Dyson*, this in- 
determinacy is avoided by observing that the normalization of y® (a) must 
not be altered by the electromagnetic interaction. The normalization condition 
for the electron field is, however, just the time-component of the current 
operator, and hence (A. 15) must vanish. This means 


f= 4 
and agrees in the e?-approximation with the result of Dyson.!® 
Department of Mechanics and Mathematical Physics, University of Lund. 


8 J, ScHwincER, Phys. Rev. 76, 790 (1949). 
* Cf. also R Karprius--N. M. Kron, Phys. Rev. 77, 542 (1950). 
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